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A new luminosity function for galaxies can be built starting from the 
product of two random variables X and Y represented by a gamma variate 
with argument 2 . The mean , the standard deviation and the distribu- 
tion function of this new distribution are computed. This new probability 
density function is assumed to describe the mass distribution of galaxies. 
Through a non linear rule of conversion from mass to luminosity a second 
new luminosity function for galaxies is derived. The test of reliability of 
these two luminosity functions was made on the Sloan Digital Sky Survey 
(SDSS) in five different bands. The Schechter function gives a better fit 
with respect to the two new luminosity functions for galaxies here derived. 



PACS numbers: 02.50.Cw Probability theory; 98.62.Ve Statistical and cor- 
relative studies of properties (luminosity and mass functions; mass-to-light ratio; 
Tully-Fisher relation, etc.) 



1. Introduction 



Given two independent non-negative random variables X and Y, their 
product XY represents an active field of research. When X and Y are Stu- 
dent's t random variables the product XY is applied in the field of finance pQ 
. When X and Y are n-Rayleigh distribution , the application can be the 
wireless propagation research [2] . In Section [2] this paper explores the 
product XY when X and Y are gamma variate with argument 2. Section [3] 
explores the connection between the Voronoi Diagrams and galaxies. Sec- 
tion H] reports two new luminosity functions for galaxies as deduced from 
the product XY. 
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2. The new distribution of probability 

The starting point is the probability density function ( in the following 
PDF) in length , s , of a segment in a random fragmentation 

p(s) = Aexp (— Xs)ds , (1) 

where A is the hazard rate of the exponential distribution. Given the fact 
that the sum , u , of two exponential distributions has PDF 

p(u) = A 2 wexp (— \u)du . (2) 

The PDF of the ID Voronoi segments , I, ( the midpoint of the sum of two 
segments) can be found from the previous formula inserting u = 21 

p{l) =2AZexp(-2AZ)cZ(2A0 . (3) 

On transforming in normalized units x = i we obtain the following PDF 

p(x) = 2xexp(-2x)d(2x) . (4) 

When this result is expressed as a gamma variate we obtain the PDF (for- 
mula (5) of [3]) 

H(x;c) = ^-^(cxf" 1 exp(-cx) , (5) 

where < x < oo , c >0 and T(c) is the gamma function with argument c; 
in the case of ID Voronoi Diagrams c = 2. It was conjectured that the area 
in 2D and the volumes in 3D of the Voronoi Diagrams may be approximated 
as the sum of two and three gamma variate with argument 2. Due to the 
fact that the sum of n independent gamma variates with shape parameter 
Cj is a gamma variate with shape parameter c = Ylf c i the area and the 
volumes are supposed to follow a gamma variate with argument 4 and 6 
[HE]. This hypothesis was later named "Kiang's conjecture", and equation 
([5]) was used as a fitting function , see [BJ [7] , or as an hypothesis to accept 
or to reject using the standard procedures of the data analysis, see [SJ [9]. 
PDF (JSJ can be generalized by introducing the dimension of the considered 
space, d(d = 1, 2, 3), 

H(x;d) = -^—(2dx) 2d - 1 exp(-2dx) . (6) 



Two other PDF are suggested for the Voronoi Diagrams: 
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• The generalized gamma function with three parameters (a, b, c) 

1)0./ c 

f(x;b,c,d) =c^j-^x a 1 exp (-6x c ) , (7) 

see [9], [EH- 

• A new analytical PDF of the type 

FN(x;d) = Const x x^r 1 exp (- (3d + l)x/2) , (8) 

where 

Const = ——^5 ics J — , (9) 

2 2 3 / 2 d (3 d + I) -372 d r (3/2 d + 1/2) 

and d(d = 1,2,3) represents the dimension of the considered space , 
see [UJ. 

Experimentally determined physical quantities are usually derived from 
combinations of measurements, each of which may be considered a ran- 
dom variable subject to a known distribution law. The distribution law of 
the sought-for physical quantity, however, is generally not known or deter- 
minable analytically, except for a linear superposition of random variables 
, i.e. the sum of n independent gamma variates, [12J. Physical quantities 
represented as products of random variables are of especial interest. Com- 
puter simulations of products of normally distributed random variables , 
as an example , lead to distributions that are not normal, in some cases 
markedly so with pronounced skewness, depending upon the parameters of 
the component distributions, [12]. Consider , for example, the product of 
two random variables X ~ N(0, 1) and the random variable Y ~ N(0, 1), 
the PDF of V = XY is , see [13] , 



h{v) = { * s '9 num ^ ~oo<v<0 

[ Kq(v * signum{v))/iT < v < oo . 

This PDF has a pole at v = . 

We now explore the product of two gamma variate with argument 2. 
We recall that if X is a random variable of the continuous type with PDF , 
f(x), which is defined and positive on the interval < x < oo and similarly 
if Y is a random variable of the continuous type with PDF g(y) which is 
defined and positive < y < oo , the PDF of V = XY is 

h(v)= [°° g( V -)f{x)-dx . (11) 
Jo X X 



4 



PRODUCT PRINTED ON JUNE 25, 2008 



Here the case of equal limits of integration will be explored , when this is 
not true difficulties arise [303] . When f(x) and g{y) are gamma variates 
with argument 2 the PDF is 

/•oo i« p -2x — 2 | 

h(v)= — dx = 32vK Uy/v) , (12) 

Jo x 

where K u (z) is the modified Bessel function of the second kind |X5|, [16] with 
v representing the order, in our case 0. The distribution function ( in the 
following DF) is 

16 v 2 K (4y^) 2F1 (l,2;3;4w) + 16v 5 / 2 K 1 (4y^) 2?i (1, 3; 3; 4 v) , 

(13) 

where 2-^1(0,6; c; v) is a regularized hypergeometric function, see [15| [T71 

EH]. 

The mean of the new PDF, h(v), as represented by formula (|12p is 

/■oo 

<v >= v x 32uiT (407) dv = 1 , (14) 
Jo 

and the variance 

/■oo _ r 

a 2 = (v- I) 2 x 32vK U^)dv = - . (15) 
Jo 4 

The mode , m , is at v =0.15067 . Figure [T] reports our function /i(x) as 
well the Kiang function H(x; c) for three values of c. Asymptotic series are 

h{v) 16 (2 In (2) + In (u) + 2 7 )?j K<1 , (16) 
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where 7 is the Euler-Mascheroni constant. 

3. Voronoi Diagrams and galaxies 

The applications of the Voronoi Diagrams, see [19] and [20] . in astro- 
physics started with [3] where through a Monte Carlo experiment, the area 
distribution in 2D and volume distribution in 3D were deduced. The ap- 
plication of the Voronoi Diagrams to the distribution of galaxies started 
with [21], where a sequential clustering process was adopted in order to in- 
sert the initial seeds. Later a general algorithm for simulating one-dimensional 
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Fig. 1. Plot of h(x) as function of x (full line) , H (x ;c ) when c — 2 (dashed), H 
(x ;c ) when c = 4 (dot-dash-dot-dash) and H (x ;c ) when c = 6 (dotted). 

lines of sight through a cellular universe was introduced [22] . The large mi- 
crowave background temperature anisotropies over angular scales up to one 
degree were calculated using a Voronoi model for large-scale structure for- 
mation in [23J and [24] . The intersections between lines that represent the 
'pencil beam' surveys and the faces of a three-dimensional Voronoi tessella- 
tion has been investigated by [25] where an exact expression is derived for 
the distribution of spacings of these intersections. Two algorithms (among 
others) that allow to detect structures from galaxy positions and magnitudes 
are briefly reviewed : 

• A Voronoi Galaxy Cluster Finder (VGCF) that uses galaxy positions 
and magnitudes to find clusters and determine their main features: 
size, richness and contrast above the background , see [26\ 127]. 

• An automated procedure for structure finding, involving the Voronoi 
tessellation allows to build a catalogue (called PF) of galaxy structures 
(groups and clusters) in an area of 5,000 square degrees in the southern 
hemisphere, see [2HJ [2j|] . 

This section first explores how the fragmentation of a 2D layer as due 
to the 2D Voronoi Diagrams can be useful or not in describing the mass 
distribution of galaxies. The large scale structures of our universe are then 
explained by the 3D Voronoi Diagrams in the second section. 
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3.1. Mass distribution 

Here we Analise the fragmentation of a 2D layer of thickness which 
is negligible with respect to the main dimension. A typical dimension of 
the layer can be found as follows. The averaged observed diameter of the 
galaxies is 

D°te « 0-6D^ s ax = 2700^ = 27 Mpc , (18) 

where = 4500 corresponds to the extension of the maximum void 

visible on the CFA2 slices. In the framework of the theory of the primordial 
explosions ,see [30] and [31], this means that the mean observed area of a 
bubble ,A°te, is 

j~\obs 

Aobs ~ 4 ^(±^)2 = 2290Mpc 2 . (19) 
The averaged area of a face of a Voronoi polyhedron , Ay, is 



j^obs 

A V == , (20) 

where Np is the averaged number of irregular faces of the Voronoi polyhe- 
dron, i.e. Np=l6 , see [32j [7]. The averaged side of a face of a irregular 
polyhedron , Ly , is 



L v » V A° bs m 12 Mpc . (21) 

The thickness of the layer , 5 , can be derived from the shock theory , 
see [33], and is 1/12 of the radius of the advancing shock , 

r~\obs 

*- max -X.YlMpc . (22) 



2 x 12 

The number of galaxies in this typical layer , Ng, can be found by multi- 
plying n* « 0.1 , the density of galaxies , by the volume of the cube of side 
12 Mpc : i.e. N G ps 172. 

A first application of the new PDF, h(v), as represented by formula (|12|) 
, can be a test on the area distribution of the Voronoi polygons , see Figure[2J 
Does the area distribution of the irregular polygons follow the sum or the 
product of two gamma variates with argument 2 1. In order to answer 
this question we fitted the sample of the area with h{v) , the new PDF , 
with a gamma variate with argument 4 and with a gamma variate with the 
argument as deduced from the sample. The results are reported in Table [TJ 
From a careful inspection of Table [T] it is possible to conclude that the area 
distribution of the irregular Voronoi polygons is better described by the sum 
of two gamma variates with argument 2 rather than by the product. 
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Table 1. The x 2 of data fit when the number of classes is 10 for three PDF 



PDF 



h(x) 

H(x; c) when c = 4 

H(x; c) when c = 3.7 

FN(x; d) Ferenc & Neda formula 



X 



62.8 
28.3 
23.55 

when d = 2 20.2 




Fig. 2. The Voronoi-Diagram in 2D when random seeds are used. The selected 
region comprises 102 seeds. 



3.2. Spatial dependence 

Our method considers a 3D lattice made of pixels^ points : present 
in this lattice are N s seeds generated according to a random process. All 
the computations are usually performed on this mathematical lattice; the 
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Fig. 3. Portion of the Voronoi-Diagram V p (2, 3) when random seeds are used; cut 
on the x-y plane . For astronomical purposes we only report a slice 130° long. 
The parameters are pixels = 800 , N — 1900 and side — 2 x 16000 Km/sec 



conversion to the physical lattice is obtained by multiplying the unit by 
5 = pixds-i > wnere s ide is the length of the square/cube expressed in the 
physical unit adopted. The tessellation in 3ft 3 is firstly analyzed through a 
planar section . Given a section of the cube (characterized , for example, by 
k = the various Vi (the volume belonging to the seed i) may or may not 
cross the little cubes belonging to the two dimensional lattice . Following 
the nomenclature introduced by [32] we call the intersection between a plane 
and the cube previously described as V p (2, 3); a typical example is shown in 
Figure [3j 

For astronomical purposes is also interesting to plot the little cubes 
belonging to a slice of 6° wide and about 130° long, see Figure HI 
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cz (km/s) 



Fig. 4. Polar plot of the little cubes belonging to a slice 130° long and 6° wide. 
Parameters as in Figure [3l 



The new PDF as given by formula (|12[) can represent a PDF in mass for 
the galaxies. The luminosity function , in the following LF, for galaxies is 
then deduced by introducing a linear and a non linear relationship between 
mass and luminosity. A special section is devoted to the parameters deter- 
mination of the two new LF for galaxies. The dependence of the number of 
galaxies with the redshift is then analyzed by adopting the first new LF. 



where a sets the slope for low values of L , L* is the characteristic luminosity 
and <£* is the normalization. This function was suggested by [33] in order to 
substitute other analytical expressions, see for example, formula (3) in |35j . 
Other interesting quantities are the mean luminosity per unit volume, j , 



4. A new luminosity function for galaxies 



4-1- Schechter luminosity function 
A model for the LF of galaxies is the Schechter function 




(23) 




(24) 



and the averaged luminosity ,(L) 



(25) 
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Table 2. The parameters of the Schechter function from 2dFGRS , SDSS , MGC 
and CFA . 



parameter 
M* [mags] 
a 

[h Mpc- 3 } 



2dFGRS 
-19.75 ±0.05 
-1.09 ±0.03 
(2.02±0.02)10~ 2 



SDSS (r*) band 
-20.83 ± 0.03 
-1.2 ±0.03 
(1.46±0.12)10~ 2 



MGC 
-19.60 ±0.04 
-1.13 ±0.02 
(1.77 ±0.15)10- 2 



CFA 
-18.79 ±0.1 
-l.±0.1 
(4.0±0.1)10~ 2 



where T is the gamma function and its appearance is explained in [36j . 
An astronomical form of equation (|23|) can be deduced by introducing the 
distribution in absolute magnitude 

$>{M)dM = (0.4ZnlO)$*10 a4(Q+1)(M *- M) 

x exp(-10 a4 ( M *- M VM , (26) 

where M* is the characteristic magnitude as derived from the data. At 
present this function is widely used and Table [2] reports the parameters 
from the following catalogs 

• The 2dF Galaxy Redshift Survey (2dFGRS) based on a preliminary 
subsample of 45000 galaxies, see [57] . 

• The r*-band LF for a sample of 11,275 galaxies from the Sloan Digital 
Sky Survey (SDSS) , see [38J. 

• The galaxy LF for a sample of 10095 galaxies from the Millennium 
Galaxy Catalogue (MGC) , see [39]. 

• The CFA Redshift Survey [ID] that covered 9063 galaxies with Zwicky 
m magnitude < 15.5 to calculate the galaxy LF over the range 13 < 
M < 22. 

Over the years many modifications have been made to the standard Schechter 
function in order to improve its fit: we report three of them. When the fit 
of the rich clusters LF is not satisfactory a two-component Schechter-like 
function is introduced , see [H] 

<£* L L 
L max > L > L Dwarf : <S>{L)dL = (—)(—)" e xp(-— )ciL , 

(27) 

L Dwarf >L>L mm : ^L)dL = (^^)(-^—r^fdL , 

Li ^Dwarf 
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where 

wart \a wart \ 

^Dwarf = $ ( ) eX P ( . 

This two-component function defined between L max and L m i n has two ad- 
ditional parameters: Lf) war f which represents the magnitude where dwarfs 
first dominate over giants and an war f the faint slope parameter for the 
dwarf population. 

Another function introduced in order to fit the case of extremely low 
luminosity galaxies is the double Schechter function , see [42] : 



<$>{L)dL = ^exp(-L/L*) 



L\ ai / L\ a2 



(28) 



where the parameters $* and a which characterize the Schechter function 
have been doubled in and 0^2- 

4-2. A linear mass-luminosity relationship 

We start by assuming that the mass of the galaxies , Ai, is distributed 
as h{M). We then assume a linear relationship between mass of galaxy and 
luminosity , L , 

L = 1ZM , (29) 

where 7Z represents the mass luminosity ratio ~ (10 — 15) , see |43| . When 
L* represents the scale of the luminosity. Equation (|12p changes to 

32LK U^=) L 

*(L)dL = **x L ; d T* ■ (30) 

where is a normalization factor which defines the overall density of galax- 
ies , a number per cubic Mpc. The mathematical range of existence is 
< L < oo. The mean luminosity per unit volume, j , 



J 



/•CO 

/ Lm(L)dl = L*^>* . (31) 
Jo 



The relationship connecting the absolute magnitude, M , of a galaxy 
with its luminosity is 



1O °- 4 (^ oi ,0- M ) , (32) 



where M bo i q is the bolometric luminosity of the sun , which according to 
is M M ' =4.74 . 
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Table 3. The full range in magnitudes , the selected range in magnitudes , the 
parameters of our function ([33]) , % 2 , AIC and BIC of our function and the 
Schechter function for the SDSS catalog 



parameter 


u* 


9* 


7* * 


i 


z 


full range [mags] 


-22,-15.8 


-23.4,-16.3 


-24.48,-16.3 


-24.5,-17.2 


-23.7,-17.^ 


selected range [mags] 


-20.65,-15.8 


-22.09,-18.2 


-22.94,-18.5 


-23.42,-18.5 


-23.7,-19 


M*[mags] 


-17.23 


-18.74 


-19.63 


-20.05 


-20.37 


vp* 


0.052 


0.033 


0.028 


0.027 


0.026 


x 2 


563 


1151 


2758 


4202 


4588 


AIC, k = 2 


567 


1155 


2762 


4206 


4592 


BIC, k = 2 


575 


1163 


2770 


4215 


4601 


X 2 Schechter 


330 


456 


1497 


1916 


2694 


AIC Schechter, k = 3 


336 


462 


1503 


1922 


2700 


BIC Schechter, k = 3 


349 


474 


1515 


1935 


2713 



A more convenient form in terms of the absolute magnitude M is 
^{M)dM = 12.8 <T io°- 8M *-°- 8M ir ( 4 . 010 o.2M*-o.2M^ ln ( 10 ) dM (33) 

This data oriented function contains the parameters M* and Vl/* which can 
be derived from the operation of fitting the observational data. 

In order to make a comparison between our LF and the Schechter LF we 
first down-loaded the data of the LF for g alaxies in the five bands of SDSS 
available at |http:/ /cosmo.nyu.edu/blanton/lf.html[ The LF for galaxies as 
obtained from the astronomical observations ranges in magnitude from a 
minimum value , M m i n , to a maximum value , M max ; details can be found 
in [15] and For our purposes we then introduced an upper limit , Mu m 
, for the absolute magnitude in order to check the range of reliability of our 
LF as represented by equation (|3"3"j) . It is interesting to stress that Mu m is 
used only for internal reasons and is connected to how the LF, as a function 
of the absolute magnitude, reaches the maximum. Table reports the 
original range in magnitude of the astronomical data , the selected range 
adopted for testing purposes , the three parameters of our function , the \ 2 
of the fit and the \ 2 OI th e Schechter function for the five bands of SDSS. 
The Schechter function , the new function and the data are reported in 
Figure [5], Figure [6], Figure [7], Figure [8], and Figure [9] when the u* ,g* ,r* , 
i* and z* bands of SDSS are considered. 
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absolute magnitude M 

Fig. 5. The LF data of SDSS (it*) are represented through the error bar. The fitting 
continuous line represents our LF (f3"3"| and the dotted line represents the Schechter 
function. 




1 -22 -21.5 -21 -20.5 -20 _ -19.5 -18 -18.5 

absolute magnitude M 

Fig. 6. The LF data of SDSS(g*) are represented through the error bar. The 
fitting continuous line represents our luminosity function (I33[) and the dotted line 
represents the Schechter function. 
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absolute magnitude M 

Fig. 7. The LF data of SDSS(r*) are represented through the error bar. The fitting 
continuous line represents our LF (f3"3"| and the dotted line represents the Schechter 
LF. 




absolute magnitude M 

Fig. 8. The LF data of SDSS(i*) are represented through the error bar. The fitting 
continuous line represents our LF ((33)) and the dotted line represents the Schechter 
LF. 
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-23 -22 -21 -2 

absolute magnitude M 



-20 



-19 



Fig. 9. The LF data of SDSS(z*) are represented through the error bar. The 
fitting continuous line represents our luminosity function (|33[) and the dotted line 
represents the Schechter LF. 



4-3. A non linear mass-luminosity relationship 

Also here we assume that the mass of the galaxies , Ai, is distributed 
as h{M). The first transformation is 



M [ ] - 

L* 



(34) 



where L is the luminosity , 1/a is an exponent that connects the mass to 
the luminosity and L* represents the scale of the luminosity. Equation (112p 
changes to 



att -2 + a 2 / 1 1 \ 

y NL (L)dL = — 32L —L*~«K Ul^L*~^) d- 



(35) 



a v / L* 

where is a normalization factor and the apex NL stands for nonlinear. 
The mean luminosity per unit volume, j , 

POO 

j= L^ NL (L)dl = A- a L*^*(T(2 + a)) 2 . (36) 
Jo 

The second transformation connects the luminosity with the absolute mag- 
nitude 



^ NL {M)dM = 12.8 V* 10 



* m-0.8 



-M*+M 



Kg [ 4.010 



0.2 ~ M *+ M \ In (10) 



-dM (37) 
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Table 4. The full range in magnitudes , the selected range in magnitudes , the 
parameters of our mass-luminosity function (|37|) , x 2 , AIC and BIC of our mass- 
luminosity function and the Schechter LF for the SDSS catalog 



parameter 


u* 


9* 




i 


z 


full range [mags] 


-22,-15.8 


-23.4,-16.3 


-24.48,-16.3 


-24.48,-17.2 


-23.7,-17.48 


selected range [mags] 


-20,-15.78 


-22,-18.2 


-22.94,-18.5 


-23.42,-19.3 


-23.7,-20 


a 


0.98 


0.95 


1.07 


1.04 


1.05 


M* [mags] 


-17.27 


-18.85 


-19.47 


-19.98 


-20.28 


\J/* 


0.05 


0.03 


0.033 


0.027 


0.027 


x 2 


552 


803 


1180 


306 


475 


AIC, k = 3 


558 


809 


1186 


312 


481 


BIC, k = 3 


570 


821 


1199 


325 


493 


X 2 Schechter 


330 


456 


1497 


1863 


2292 


AIC Schechter, k = 3 


336 


462 


1503 


1869 


2298 


BIC Schechter, k = 3 


349 


474 


1515 


1882 


2310 



Table 5. The parameters of the first LF 

based on data from CFA Redshift Survey ( triplets generated by the author) 



CFA 


M*[mags] 


-19 ±0.1 


^>*[h Mpc~ 3 ] 


0.4 ±0.01 



The parameters that should be deduced from the data are M* , a and 
Table H] reports the original range in magnitude of the astronomical data , 
the selected range adopted for testing purposes , the three parameter of our 
function , the \ 2 °f the fit and the \ 2 °f the of the Schechter function for the 
five bands of SDSS. Also here in the u* case the astronomical range and the 
selected range are coincident. In the absence of observational data which 
represent the LF , we can generate them through Schechter's parameters, 
see Table [2j this is done, for example, for the CFA Redshift Survey ,see [4*0] . 
The parameters of the first LF ( equation (|33p ) are reported in Table [5] 
where the requested errors on the values of luminosity are the same as the 
considered value. 
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4-4- Parameters determination 

The theoretical LF for galaxies can be represented by an analytical func- 
tion of the type <&(M*, <j>*,Pz) where M*,(p* and p^ represent the scaling 
magnitude , the number of galaxies per unit Mpc and a generic third pa- 
rameter. Once the observational data are provided in n triplets made by 
absolute magnitude , 4> as tr ( in units of number per h~~ 3 Mpc~ 3 per mag) and 
00 ( the error on ) we can deduce these three parameters in the following 
ways. 

• A scanning on the presumed values of the parameters that are un- 
known. The three parameters are those that minimize the merit func- 
tion x 2 computed as 

X 2 = E(^^) 2 • (38) 



• A nonlinear fit through the Levenberg-Marquardt method ( subrou- 
tine MRQMIN in [TB] ) . In this case the first derivative of the LF with 
respect to the unknown parameters should be provided. 

Particular attention should be paid to the number of parameters that are 
unknown : two for the new LF as represented by formula (|33p , three for 
the Schechter function (formula (|26p ) and the new mass-LF relationship 
(formula (|37p ). The Akaike information criterion {AIC) , see [37] , is defined 

as 

AIC = 2k- 2ln(L) , (39) 

where L is the likelihood function and k the number of free parameters of the 
model. We assume a Gaussian distribution for the errors and the likelihood 
function can be derived from the x 2 statistic L oc exp(— ^-) where x 2 has 
been computed trough equation ([38]) . see [3E],[39]. Now AIC becomes 

AIC = 2k + X 2 ■ (40) 

The Bayesian information criterion (BIC), see [SD], is 

BIC = k ln(n) - 2ln(L) , (41) 

where L is the likelihood function , k the number of free parameters of the 
model and n the number of observations. 
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4-5. Dependence from the redshift 

The joint distribution in z (redshift) and / (flux) for galaxies , see for- 
mula (1.104) in [51] or formula (1.117) in [43] , is 

dN , c , k a T , z 2 



dfldzdf v iJ v z crit 



4vr (_)V*(— ) , (42) 



where df2 , dz and <if represent the differential of the solid angle , the redshift 
and the flux respectively. The critical value of z , z cr u , is 



tt2 t * 

"crit » r 9 

4vr/cf 



(43) 



The number of galaxies , N s (z, f m ini fmax) comprised between a minimum 
value of flux, f m i m and maximum value of flux f max > can be computed 
through the following integral 

Ns(z) = / 47r(-^)V*(^)# . (44) 

Jfmin -"0 Z crit 

This integral does not have an analytical solution and we performed a nu- 
merical integration. The number of galaxies in z and / as given by for- 
mula (@2D has a maximum at z — Zp OS — max , where 

Zpos—max 1.3798 X Z cr n , (^^) 

that can be re-expressed as 



0.3892Vl0 a4M O-°- 4M X 

VfcL ' ( } 

where Mq is the reference magnitude of the sun at the considered bandpass, 
Hq is the Hubble constant and cl is the velocity of the light. 

From the point of view of the astronomical observations the second CFA2 
redshift Survey , started in 1984, produced slices showing that the spatial 
distribution of galaxies is not random but distributed on filaments that 
represent the 2D projection of 3D bubbles. We recall that a slice comprises 
all the galaxies with magnitude m& < 16.5 in a strip of 6° wide and 
about 130° long. One such slice (the so called first CFA strip) is visible 



at the following address http://cfa- www. harvard. edu~huchra/zcat/ and is 



reported in Figure [10]; more details can be found in [52J. 

Figure [11] reports the number of observed galaxies of the second CFA2 
redshift catalogue for a given magnitude and the theoretical curve as repre- 
sented by formula ([4*2]) . 
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Fig. 10. Polar plot of the real galaxies (green points) belonging to the second CFA2 
redshift catalogue. 

The total number of galaxies of the second CFA2 redshift catalogue is 
reported in Figure [12] as well as the theoretical curve as represented by the 
numerical integration of formula (|44p . 

A typical polar plot is reported in Figure [T51 once the number of galaxies 
as a function of z is computed through the numerical integration of for- 
mula (|44|) ; it should be compared with the observations , see Figure [TUJ 



5. Summary 

The PDF of the product of two independent random variables X and Y 
as represented by two gamma variates with argument 2 has been analytically 
derived. 

The mean , the variance and the DF of this new PDF are computed. As 
an application we assumed that the mass of the galaxies behaves in the same 
way as this new PDF. The LF for galaxies can therefore derived assuming 
a linear or nonlinear relationship between mass and luminosity: in the first 
case we have a two parameter LF and in the second case a three parameter 
LF; recall that the Schechter LF for galaxies has three parameters. The 
parameter a that characterizes the non-linear relationship between mass 
and luminosity of the second LF , see equation (|37p , is found to be around 
1. 

The comparison between the two LF for galaxies here derived is per- 
formed on the SDSS data and can be done only by introducing an upper 
limit in magnitude , Mn m , in the five bands analyzed. The three tests of re- 
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Fig. 11. The galaxies of the second CFA2 redshift catalogue with 15.08 < mag < 
15.27 or 48776 M ® 2 < / < 58016 M ® 2 ( with mag representing the relative magni- 
tude used in object selection), are isolated in order to represent a chosen value of 
m and then organized in frequencies versus heliocentric redshift , (empty circles); 
the error bar is given by the square root of the frequency. The maximum in the 
frequencies of observed galaxies is at z = 0.02. The theoretical curve generated 
by the z-dependence in the number of galaxies (formula (|42[) and parameters as in 
column CFA of Table [5]) is drawn (full line) . 



liability here adopted show that the Schechter function always has a smaller 
X 2 , AIC and BIC with respect to the two new LF for galaxies here derived 
, see Tableland Tabled 

The theoretical number of galaxies as a function of the red-shift presents 
a maximum that is a function of a and / for the Schechter function; con- 
versely, when the first new LF here derived is considered , the maximum is 
a function only of /, see equation (|45p . The first new LF for galaxies , once 
implemented on a 3D Voronoi slice , allow us to reproduce the large scale 
structures of our universe. 
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